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Self-consistent approach for thermodynamics of a simplified 

pseudospin-electron model 
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We present the method of the self-consistent calculation 
of thermodynamical and correlation functions. This approach 
is based on the GRPA (generalized random phase approxima- 
tion) scheme with the inclusion of the mean field corrections. 
Investigation of a pseudospin-electron model (PEM) within 
the framework of the presented method shows that interac- 
tion between the electron and pseudospin subsystems leads 
to the possibility of either first or second order phase transi- 
tions between different uniform phases (bistability) as well as 
between the uniform and the chess-board ones. In the regime 
n — const, an instability with respect to phase separation in 
the electron and pseudospin subsystems can take place. 

Key words: pseudospin-electron model, chess-board phase, 
phase transitions, local anharmonicity, high-T c superconduc- 
tors 

PACS numbers: 71.10.Fd, 71.38. +i, 77.80.Bh, 63.20.Ry 



1 Introduction 

The crystals with the high-temperature superconductiv- 
ity (HTSC) investigated intensively during the last ten 
years possess a wide spectrum of the interesting physi- 
cal properties. The variety of effects which are realized 
only separately in the case of other types of crystals is 
characteristic of these systems. The strong electron cor- 
relation, related to the interaction of the Hubbard type 
in the conducting bands formed mainly by the supercon- 
ducting CU2-O2 plains, can be pointed out as one of the 
reasons of this unique situation. 

Another significant feature connected with the dy- 
namics of HTSC crystals is the presence of strongly an- 
harmonic elements of the structure. As it is known, it 
can be a source of instabilities of various type. 

Among the most frequently investigated HTSC crys- 
tals are the crystals of the YP^CuaCv^a group. The 
unit cell of these compounds contains, besides two su- 
perconducting planes, the chain (at the 5 <C 1 com- 
position) elements Cui-Oi, connected by CU1-O4-CU2 
bridges with conducting plains through the apical oxy- 
gen ions O4. The vibrations of these ions along the c-axis 
(perpendicularly to the plains) exhibit a strong anhar- 
monicity. Much evidence exists in support of this con- 
cept. One can mention the results of experimental in- 
vestigations (EXAFS data [1,2], Raman scattering [3-6], 



dielectric measurements [7,8]) where the two equilib- 
rium positions of O4 ion (two different values of the 
-R04-CU2 distance) have been observed that can point 
out to the presence of the local double-minimum poten- 
tial well. Similar conclusions were made based on con- 
sideration the local polaron phenomena [9,10], electron 
transfer processes through O4 ions [11,12] or bistabilitics 
in the normal phase temperature region [13]. 

Besides, a connection between positions of O4 ions 
and the energy of electron states in CU2-O2 plains plays 
an important role in YBa2Cu307_a crystals. The data 
given in [14] point to the existence of a significant cor- 
relation between the occupancy of electron states of the 
Cu2 ion and the Ro 4 -Cu 2 distance as well as to the de- 
crease of this distance at the transition from the metal- 
lic orthorombic phase to the semiconducting phase (that 
takes place at S > S* = 0.55). These and other similar 
facts suggest the presence of a large electron- vibrational 
coupling. 

By now, the description of the locally anharmonic 
subsystem in HTSC crystals develops with the use of two 
different approaches. The first one, which is chronolog- 
ically older, is based on the phonon anharmonic model 
(p 4 [15], or, more recently, on the model ip 3 + ip [13]. 
This approach was used while consideration the polaron 
effect [16] and also while investigating the effect of anhar- 
monicity on the superconducting transition temperature 
and on the possibility of the modulated (of the CDW 
type) phase creation [17]. 

In the second approach, which is more appropriate at 
the strong electron-vibrational coupling and in the cases 
when two equilibrium positions of anharmonic ions re- 
ally exist, the pseudospin formalism is used; the pseu- 
dospin variable Sf = ±1/2 defines these two positions. 
This scheme being applied to the YBa2Cu3 07_<5 type 
systems started from works [18] and [19]. 

On the basis of the second model which was called 
as pseudospin-electron model (PEM) a possible connec- 
tion between the superconductivity and lattice instabil- 
ity of the ferroelectric type in HTSC has been discussed 
[19,20]. The description of the electron spectrum and the 
electron statistics of the PEM was given in [21] in the 
framework of the temperature two-time Green function 
method in the Hubbard-I approximation. 
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A series of works has been carried out in which the 
pseudospin (S Z S Z ), mixed {S z n) and charge (rm) cor- 
relation functions were calculated. It has been shown 
with the use of the generalized random phase approx- 
imation (GRPA) [22] in the limit of infinite single-site 
electron correlations (U — > oo) [23,24], that there exists 
a possibility of divergences of these functions at some 
values of temperature. This effect was interpreted as 
a manifestation of dielectric instability or ferroelectric 
type anomaly. The tendency to the spatially modulated 
charge and pseudospin ordering at the certain model pa- 
rameter values was found out. 

On the other hand, the case of absence of the term 
describing electron transfer in Cu 2 -0 2 plains (t^ = 0) 
with the inclusion of the direct interaction between pseu- 
dospins was considered within the mean field approxima- 
tion [25,26] (see also the short review in [27]). The first 
or second order phase transitions with the jumps of {S z } 
and electron concentration n values in the (i = const 
regime were obtained. An instability with respect to 
phase separation in the electron and pseudospin subsys- 
tems can take place in regime n — const. 

The analysis of ferroelectric type instabilities in the 
two-sublattice model of the apex oxygen subsystem in 
high temperature superconducting systems has been 
made [28] . The influence of oxygen nonstoichiometry on 
localization of apex oxygen in YBa2Cu30y_5 type crys- 
tals was studied in the work [29]. 

In the present work we propose the self-consistent 
scheme for calculation of mean values of pseudospin and 
electron number operators, grand canonical potential as 
well as correlation functions for the case of the U = 
limit (the simplified PEM). The approach is based on 
the GRPA with the inclusion of the mean field type con- 
tributions coming from effective pseudospin interactions 
via conducting electrons [30] . The main attention is paid 
to the thermodynamics of phase transitions. The possi- 
bilities of phase separation and chess-board phase ap- 
pearance are investigated. 

2 Pseudospin-electron model 

Hamiltonian of the pseudospin-electron model has the 
following form: 

H=H + ^ kjbtrbjv, 

ija 

H = ^2 \ Un^rm -fi n i<T +g n irT Sf -hS- > , 

i \ a a ) 

where the strong single-site electron correlation U, in- 
teraction with the anharmonic mode (c/-term) and the 



energy of the anharmonic potential asymmetry (h-tevm) 
are included in the single-site part; [i is the chemical 
potential. The second term in the Hamiltonian describes 
the electron hopping from site to site (the electron trans- 
fer parameter tij). 

The formalism of electron annihilation (creation) op- 
erators a ia = b lrT P+, a la = b ia P~ (P^ = \ ± Sf) act- 
ing at a site with the certain pseudospin orientation is 
introduced. The calculation is performed in the strong 
coupling case (g » t) using of single-site states as the 
basic one. 

Ha=^2{e(n l ^+n ll ) - hSf}, (1) 

i 

#int= t v( a ta a ]°+ a i<J h ]°+ h t<j a ]<?+ h t<j h ]<?)- 
ija 

Here e = —fi + g/2, e = — fj, — g/2 are energies 
of the single-site states. We consider here the simplified 
PEM (the case U = 0). 

Expansion of the calculated quantities in terms of 
electron transfer leads to the infinite series of terms con- 
taining the averages of the T-products of the , cn a op- 
erators. The evaluation of such averages is made using 
the corresponding Wick's theorem. The results are ex- 
pressed in terms of the products of nonperturbed Green 
functions and averages of a certain number of the pro- 
jection operators Pf~ which are calculated by means of 
the semi-invariant expansion [30]. 

Single-electron Green function (calculated in Hub- 
bard-I type approximation) may be written as the fol- 
lowing chain diagram: 




(2) 



In the adopted approximation the diagrammatic se- 
ries for the pseudospin mean value can be presented in 
the form 





(3) 



Here we use the following diagrammatic notations: 

P+ P.~ 



0= S ii — C=€>- =9i(un) = - 1 h : 1 — =, nonper- 

turbated electron Green function — — ={gi(u n )) , 
wavy line is the intersite hopping tij . Semi-invariants are 
represented by ovals and contain the <5-symbols on site in- 
dexes. In the spirit of the traditional mean field approach 
[30] the renormalization of the basic semi-invariant by 
the insertion of independent loop fragments is taken into 
account in (||). 
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The diagrammatic series for the electron concentra- 
tion mean value is the following: 




(4) 



1 



(P?) 



where I I =rij 

lW n — fc~ lLU n - fc- 

e<*=(e,i),P?=(Pt,Pr). 

In the same approximation the grand canonical po- 
tential and pseudospin correlation function are, respec- 
tively: 



Ail = 




(5) 



(6) 



First term in equation (^) takes into account a direct 
influence of the internal effective self-consistent field on 
pseudospins: 



(7) 



leading to the renormalization of the second-order semi- 
invariant due to the inclusion of "single-tail" parts. Sec- 
ond term in equation (^|) describes an interaction be- 
tween pseudospins which is mediated by electron hop- 
ping. 




Hence, with respect to the GRPA scheme an influ- 
ence of the internal effective self-consistent field on pseu- 
dospins is taken into account by means of the inclusion of 
the mean field type contributions into the expressions for 
all thermodynamic quantities. In the approach presented 
the correlation functions are calculated consistently with 
thermodynamical functions. 

The consistency of the expressions (U)-(U) can be 
checked explicitly [30] using the thermodynamical rela- 
tions: 



dfl 



= (S z 



d(ph) 



{S z S z ) q=0 



At high temperatures eq. (||) possesses only the uni- 
form solution (Sf) = (S z ). The phase transitions be- 
tween uniform phases with different pseudospin mean 
values (S z ) were analyzed in [30]. For the first time the 
possibility of the dielectric instabilities has been done 
for the PEM in the limit of the strong electron corre- 
lation (U -> oo) in [23,24]. Results of ref. [30], where 
the opposite case of U — > was analyzed, are in good 
agreement with the exact ones for the U — > PEM in the 
limit of infinite spatial dimension [31]. At the same time, 
a complete description of such transitions was obtained 
in [25,26] for the PEM with direct interaction between 
pseudospins (ty = limit). 

On the other hand, the solution of eq. (||) for pseu- 
dospin correlator has the form 



{S z S z ) q 



1/4- (S z ) 2 



n 



,(i-<^> 2 )' 



(8) 



where II q characterize an interaction between pseu- 
dospins via electron subsystem: 



n 



(9) 



and its singularities 



n 



(i 

<? v 4 



z\2\ 



(s z ) 







gives the points of the instability of uniform phase with 
respect to fluctuations with wave vector q. The typi- 
cal dependence of the temperature of instability on the 
fluctuation wave vector q for square lattice is shown in 
figure 1 and one can see that for some model parameter 
values uniform phase become unstable with respect to 
fluctuations with q — (tt, it) (chess-board phase). 
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Figure 1. Temperature of the uniform phase instability ver- 
sus the fluctuation wave vector q = (q, q) for square lattice 
(h = 0.272, (j, = -0.36, W = 0.2, g = 1, the lattice spacing is 
seted to 1): W is the halfwidth of the initial electron band. 



So, below we shall take into account the possibility of 
the chess-board phase appearance in our analysis. 

3 The chess-board phase 

To take into account the modulation of the pseudospin 
and electron distribution we introduce two kinds of sites: 
the (Sf) corresponds to the one sublattice, and the (S 2 ) 
to the other one. The nearest-neighbour hopping exists 
only between the sublattices. 

Single-electron Green functions (||) in this case are 
equal: 



G\ k {io n ) 

G2fc(w„) 



1 - t 2 k gi{ijj n )g 2 {u>r, 



(10) 



52 (w„) 



1 - t 2 gi(uj n )g 2 (u n ) 



where gi(uj n ) and g 2 (w„) are nonperturbated Green 
functions for the sublattices 1 and 2 respectively. The 
poles of functions G k {u> n ) determine the single-electron 
spectrum. The equation for the spectrum has the form: 



x 4 - (g 2 /2 + t\)x 2 - gtU(Sl) + (S^))x + 
+g 4 /16-g 2 t 2 (S z 1 ){S^=0, 



(11) 



where x = iuj n + /i. 

The roots £i(i fc ) ^ e 2 (t k ) ^ e 3 (t k ) ^ s 4 (t k ) of the 
equation ( |Tl| ) form four subbands. The widths of sub- 
bands depend on the mean values of pseudospins. 

The branches £\(t k ), £ 2 (t k ) on the one side and 
£3(^)7 £i{tk) on the other one coincide at t k = 



{si, 2 (t k = 0) = 3/2, e 3A {t k = 0) = -g/2) and form 
two pairs of bands which are always separated by gap. 

The equations for pseudospin mean values (^) can be 
written now in the form: 



1 ( g i o_ e -/3e 

(S?) = - tanh )P(h+a\-a\) + In 



(12) 



I = 1,2; where expressions for the effective self-consis- 
tent fields are 



i= jv E **( £ - e ~) E (**)-/*]> 



(13) 



i=l 



A 



ei(t k ) + g{Sf,) 



{£i(tk)-ej(tk))(ei(tk)-e p (tk))(ei(th)-£ m (th)) , 
Expression for electron number mean value follows from 



i+n 2 ) = —^2^2n[ei{t k )-fi}- 



(14) 



k i=l 



((P 1 + )+(P 2 + ))n(f) + ((Pr) + (P 2 -)Hs) 



and grand canonical potential (||) can be written in the 
analytical form, also: 



An 



N(3 

E 



— Vln 



II COSh 



i=l 



,. (cosh f e)2 (cosh §£) 2 



+ 



(15) 



(=1,2 



E 

1=1,2 



1 f /3 

--lncosh|-(/i+4-Q!i) 



In 



l+e"^ 
1+e-^ 



- — In cosh < — /i 
P I 2 



In 



1+e 



-0s 



4 Numerical results 

In the investigation of equilibrium conditions we shall 
separate two different regimes in which the system can 
exist: 

a) the [i — const regime; it is supposed that the elec- 
tron states of other structure elements, which are not 
included explicitly into the PEM, play a role of a ther- 
mostat, that ensures a constant value of the chemical po- 
tential fi (despite the possible changes of temperature, 
field h and other characteristics of the model). In this 
case the minimum of the grand canonical potential fl is 
a condition of thermodynamical equilibrium; 
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b) the regime n = const; this situation is more cus- 
tomary at the consideration of electron systems and it 
means that the chemical potential is now the function 
of T, h etc. and depends on the electron concentration. 
The minimum of the free energy F = CI + /iiV is the 
equilibrium condition in this case. 



4.1 fi = const regime 

For this regime the equilibrium is defined by the min- 
imum condition of fl ( |15|) t hat form the equations for 
pseudospin mean values([l2|) and expression for electron 
concentration (|lj). The calculated field dependences of 
(Si — S2) (order parameter for chess-board phase) and 
grand canonical potential (which are determined by the 
solutions of the eqs. fll^)) are presented in figure 2 for 
g 3> W and low temperature. 



(.S t HS 2 ) 
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-0.32 



0.25 



0.30 




Figure 2. Field dependence of the grand canonical poten- 
tial and order parameter (T = 0.005, /i = -0.36, W = 0.2, 
g — 1). Dotted and solid lines correspond to the uniform and 
chess-board phases, respectively. 



Comparison of the grand canonical potential fl val- 
ues for uniform and chess-board phases leads to the con- 
clusion that the modulated phase is thermodynamically 
stable at intermediate values of h parameter in the re- 
gion between points a and b. These points correspond to 
the first and second order phase transition, respectively. 
The rapid jump-like change of order parameters is ac- 
companied by the rapid changes of the subbands widths 
and, as a result, electron concentration. 




Figure 3. /x - h phase diagram (T = 0.005, W = 0.2, g = 1). 
I - uniform phase, II - chess-board phase. Dashed lines - 
first order phase transitions between the uniform phases with 
different pseudospin mean values. Dotted lines - first or- 
der phase transitions between the uniform and chess-board 
phases. Solid lines - second order phase transitions. 

The transitions between uniform phases with differ- 
ent pseudospin mean values as well as corresponding 
electron concentrations (bistability), which is of the first 
order, takes place when the chemical potential \i is placed 
within the e%, £2 and partially £3, £4 bands. The tran- 
sitions between the uniform and modulated phases are 
of the first or second order and can be realized in the 
case when /1 is placed in £2 and £3 bands or between 
them. The chess-board phase exists as intermediate one 
between the uniform phases with different (S z ) and n 
values. 

The T — h phase diagram is shown in figure 4. 




The resulting phase diagram /i 
perature is shown in figure 3. 



h at the low tern- 



Figure 4. Phase diagram T — h (fJ, = -0.36, W = 0.2, 5 = 1). 
I - uniform phase, II - chess-board phase. Dashed lines - first 
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order phase transitions between the different uniform phases 
(bistability) . Dotted line - first order phase transitions be- 
tween the uniform and chess-board phases. Solid lines - sec- 
ond order phase transitions. 

With the temperature increase the first order phase tran- 
sition between the uniform and chess-board phases trans- 
forms into the first order phase transition between uni- 
form phases and, finally, disappear in the critical point 

0c- 

The diagram T — h shows the possibility of the first 
order phase transitions between uniform phases and ei- 
ther first or second order ones between the uniform and 
chess-board phases at the change of temperature. 

In figures 2 and 4 the case when the chemical poten- 
tial is placed in the lower band is presented. If the chem- 
ical potential is placed in the upper band our results 
transform according to the symmetry of the Hamilto- 
nian: /i — > — fj,, h — > 2g — h, n — > 2 — n, S z — > — S z . 



4.2 



n 



const regime 



In the regime of a fixed value of electron concentration 
the equilibrium is defined by the minimum of free energy 
F = CI + fiN. This condition form a set of equations ( |f2| ) 
and ( [l4| ) for the pseudospin mean values and chemical 
potential. The obtained dependences of F and /i on the 
electron concentration are presented in figure 5. 
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Figure 5. Dependence of the chemical potential \i on the elec- 
tron concentration n and deviation of the free energy from 
linear dependence (T = 0.005, h = 0.28, W = 0.2, g = 1). 
Dashed line - phase separation area. Dotted line - uniform 
phase. Solid line - chess-board phase. 

One can see the regions with d/i/dn ^ where states 
with a homogeneous distribution of particles are unsta- 
ble, which corresponds to the phase separation into the 
regions with different phases (the uniform and chess- 
board ones in this case) and with different electron con- 
centrations and pseudospin mean values (that is, differ- 
ent occupancies of particle positions in the anharmonic 



potential wells). In the phase separated region the free 
energy, as a function of n, deflects up and electron con- 
centrations in the separated phases are determined by 
the tangent line touch points. 

On the base of the obtained results the phase diagram 
n — h was constructed (figure 6). 




Figure 6. Phase diagram n-h(T = 0.005, W = 0.2, g = 1). 
I - uniform phase, II - chess-board phase, PS - phase sepa- 
ration area. 

The phase separation into the regions with uniform 
and chess-board phases takes place in the case when 
chemical potential is placed within the subbands £2, £3, 
that agree with the results obtained in the /1 = const 
case when within this area we had the first order phase 
transition between the corresponding phases. 

The phase separated and chess-board phase regions 
narrows with the temperature increase, but thick solid 
lines in figure 6 approach one to another faster and, for 
some temperatures, we have only the phase separation 
into the regions with uniform phases. 



5 Conclusions 

On the basis of the presented above self-consistent 
scheme for calculation of the correlation and thermo- 
dynamical functions, the energy spectrum, thermody- 
namics of phase transitions, possibility of phase sepa- 
rations as well as appearance of the chess-board phase 
have been investigated. The corresponding phase dia- 
grams were build. Such comprehensive analysis of the 
thermodynamics of the considered simplified PEM be- 
came possible due to the generalization of the traditional 
GRPA approach (in which, on the basis of calculation 
of the (S Z S Z ) correlation function, only the assumption 
about the chess-board phase appearance may be done). 
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The obtained phase diagrams remind the situation 
know for the Falicov-Kimball (FK) model (this model is 
close to PEM but differ in the thermodynamic equilib- 
rium conditions) with a rich phase diagram. However, 
contrary to this model, an existence of the phase tran- 
sitions between uniform phases is possible in our case. 
This results from the another regime of thermodynamic 
averaging (fixation of h field which is an analogous to 
the chemical potential for ions in the FK model). 

The study of the thermodynamics of the PEM within 
the framework of the presented above approach provides 
reason enough to conclude that on the basis of this model 
one can describe the phase transitions and instabilities 
in the HTSC of the YBa 2 Cu307„,5 type. Among them 
we should list the bistability effects (at the change of 
temperature), separation into phases with different elec- 
tron concentrations, pseudospin orientations (that corre- 
sponds to different localization of particles in the anhar- 
monic potential wells). In general, it is in agreement with 
the picture observed experimentally (some publications 
in this field were quoted above). 
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